Lie Algebras and Suppression of Decoherence in Open Quantum Systems 
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Since there are many examples in which no decoherence-free subsystems exist (among them all 
cases where the error generators act irreducibly on the system Hilbert space) , it is of interest to search 
for novel mechanisms which suppress decoherence in these more general cases. Drawing on recent 
work (quant-ph/0502153l we present three results which indicate decoherence suppression without 
the need for noiseless subsystems. There is a certain trade-off; our results do not necessarily apply to 
an arbitrary initial density matrix, or for completely generic noise parameters. On the other hand, 
our computational methods are novel and the result — suppression of decoherence in the error-algebra 
approach without noiseless subsystems — is an interesting new direction. 
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I. INTRODUCTION 

The central obstacle in the experimental realization of 
quantum computers has proven to be maintaining the 
quantum coherence of states Q . The main cause of this 
degradation of the quantum coherence is the coupling of 
the computer to the environment, and the decoherence 
induced by this coupling. 

A variety of schemes for protecting quantum infor- 
mation have been developed, including quantum error 
correction codes 0: S Pj i decoherence free subspaces 
0, 0, H, 1^ IToll - noiseless subsystems jllj, bang-bang 
decoupling [iJl, and topological quantum computation 
[Tsl . The first four of these techniques are closely related 
to each other and can be described in a simple unified 
framework based on representations of the algebra of er- 
rors 0, Q . More recently, Zanardi and Lloyd [l^ 
showed that topological quantum computation also falls 
into the error-algebra framework. 

Our approach, like those mentioned above, uses the 
(Lie algebraic) structure of the physical decoherence pro- 
cess itself to protect quantum coherence. We start from 
the assumption that the Lindblad operators generate a 
representation of a Lie algebra, and show that if this rep- 
resentation is irreducible, then decoherence-free subsys- 
tems and subspaces will not exist. Thus if one is to have 
any hope of controlling decoherence in the irreducible 
case, some other mechanism is needed. 



A. Markovian Dynamics 

The dynamics of a quantum system A coupled to a 
heat bath B, which together form a closed system, is 
described by a Hamiltonian : 

(iJ^^ls) + {1a® Hb) + Hi, 

where Ha, Hb, and Hj are the system, bath, and interac- 
tion Hamiltonians respectively. The Markovian Master 
equation jisl [l9l | (also "semigroup master equation" or 



SME) provides the most general form for time evolution 
of the system density operator p — pA, which acts on a 
Hilbert space Ha'- 



dp _ i 



P] + C{p) . 



(1) 



Here C is called the Lindbladian or dissipator, and takes 
the general form 



^P) 



a, 13 



Lf^,fM = [Fo.,pFl] + [F^p,Fl]. 



(2) 



(3) 



The equation of motion remains invariant under an ar- 
bitrary unitary transformation of the Lindblad opera- 
tors. The special case corresponding to C{p) =0 is the 
von Neumann equation, and describes unitary evolution. 
Therefore in general, C{p) contains precisely those terms 
in the evolution equation which are responsible for deco- 
herence. The derivation of eqns. (0) and Q from funda- 
mental assumptions is done in many places; Bacon et al. 
|20| have a particularly nice treatment. 

The coefficient matrix (oq^^) is assumed to be time- 
independent and Hermitian. By diagonalizing a, it fol- 
lows that may also be written in "diagonal standard 
form," replacing the Fa operators by suitable linear com- 
binations Va- This yields 



Xp)- 



(4) 



The noise coefficients a^^j have been absorbed within 
the operators Va- From using the relation 

La, A — [A, [p, A]] which holds for a Hermitian operator 
A, we find 

a a ^ ^ 
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If J2a — '^1 where c G M is a constant, then we have with Kraus operators 



The first term is reminiscent of a quantum channel with 
Hermitian Kraus operators. Indeed, an equivalent form 
of the Markov master equation is 
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where $ is a completely positive map 0, p. 18,11.5.5]. 

In applications to real-world systems, the matrix ele- 
ments ttap contain physical parameters such as lifetimes, 
longitudinal or transverse relaxation times, stationary 
magnetization, etc. We will show that particularly sim- 
ple dynamics emerge if the coefficient matrix is not only 
Hermitian but also real (hence symmetric) or approxi- 
mately so. 

The Fa describe various decoherence processes, and for 
this reason, they are often called error generators. The 
Fa are determined implicitly by the interaction Hamilto- 



where {Ba} is a collection of operators on Hb, called 
the heat bath operators. Let g[(7i^) denote the vector 
space of all linear transformations on Ha , not necessarily 
invertible. This is the Lie algebra of the group GL{n) 
where n = dim(7i^). We now state the main assumption 
of the paper. 

Assumption 1. There exists a Lie algebra q and a faith- 
ful representation : g — > g[(7i^) such that 



Fa^Hfa), a^Q...M 
for some linearly independent set {fa} C 0. 



(5) 



This assumption is also basic to many other studies of 
decoherence; see 0, |23, for a few. The trace of any 
generator of any representation of a compact simple Lie 
algebra is zero |23l Thm. 8.9]. Further, we are primarily 
interested in representations for which all the generators 
can be chosen to be Hermitian. This is always possible 
for compact Lie algebras [23L Sec. 2.4]. We may therefore 
take Fa to be traceless and Hermitian. 

Decoherence processes in which the Fa are Hermitian 
are mathematically much simpler than the most general 
process, and this is often the case of interest to applica- 
tions in quantum computing. We recall for clarity the 
standard example of a quantum computer made of K 
qubits with n = 2^ dimensional register Hilbert space. 

Each qubit has four possibilities: no error, or an error 
generated by one of the three Pauli matrices. This means 
that each qubit independently undergoes the action of the 
standard depolarizing channel 



P 



(6) 



Mo = v/T^ 1, M, - 



(* = 1...3). (7) 



where cr^ are Pauli matrices and Z is a normalization 
constant equal to 3/4 for the spin-1/2 representation. 

The maximum possible complexity of error generation 
is when combined errors from any number of qubits are 
generated. The 4^ — 1 error generators Fa for such a 
process are basis elements of g = 5u(n) — 5u(2^') in the 
defining (n-dimensional) representation. As the defining 
representation is irreducible, there are no DF subspaces. 
We have written (|7I) in a notation compatible with the 
generalization to arbitrary Lie algebras done in Section 
II VI and more completely in [2^ . 



B. Irreducible Representations Have No 
Decoherence-EVee Subsystems 

Knill, Laflamme, and Viola discovered a method 
for decoherence-free coding into subsystems instead of 
into subspaces which has since received much attention; 
see 1^ ■ In this section, we remark that in the 

situation of Assumption ^ existence of a DF subsystem 
implies that the representation is reducible. 

The definition of a DF subsystem begins with the ob- 
vious statement that if the Fa are Hermitian and one of 
them is the identity, then the associative algebra A they 
generate is unital and closed under adjoints. Therefore, 
C*-algebra methods may be applied. 

Let A4 {d, C) denote the space oi dx d complex matri- 
ces. Any basis for Ti. determines a matrix representation 
(j) : A M{d, C), where d = dim(7i), simply by express- 
ing each operator in this basis. The idea is that (f> may 
be repetitive^ i.e. it may happen that 
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where </)d; : A^ A4{di,C). The rii, 712, . . . label the sizes 
of the corresponding sub-matrices. In fancier notation, 
(j){a) — J2i Irii ^ 4>di{o)- This is also called the central 
decomposition. 

If this happens, then we may take a vector ^I* e 
which is also repetitive, so that the first di components 
of 'J take the form aw, where a € C, u G , and the sec- 
ond di components take the form (iv for the same v, etc. 
We may do this ni times. Then (/)(a)^' = {aw,l3w, . . .) 
where w — (f)di{a)v. The coefficient vector (a,/?, . . .) G 
can just as well be extracted after applying the op- 
erator 0(a) as before. In other words, the information 
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contained in the a, (3, ... is protected under this decoher- 
ence process. 

Theorem 1. Decoherence-free subsystems do not exist if 
the error operators Fa generate an irreducible represen- 
tation. 

Proof. If the Fa generate an irreducible matrix represen- 
tation of a semisimple Lie algebra g, then A will be the 
full algebra of Hermitian operators over H. Therefore, 
the homomorphism cj) : A ^ M{d,C), where d = dim(H) 
is not repetitive. In the above notation, ni = 1, there 
is no 712 and di = dim(7i). Second proof: noiseless de- 
grees of freedom are associated with observables in A', 
the commutant algebra. By Schur's Lemma, A' is the 
trivial algebra if (j) is irreducible. □ 

This theorem implies the corresponding result for DF 
subspaces as a special case. Conditions under which DF 
subspaces can exist have previously been studied in both 
the Lindblad (Markovian) formulation [1,0] and for the 
non-Mar kovian case Q. Theorem ^ extends their results 
to noiseless subsystems. 

C. Irreducible Representations Return to 
Equilibrium 

In this section we prove a property of the time- 
evolution of open quantum systems defined by irreducible 
representations. This property is not used directly in the 
rest of the paper, but it holds independent interest. The 
thermal Gibbs state 

is a stationary state for the Markovian dynamics of a 
system coupled to a heat bath. The question then arises: 
under what conditions does the system return to equilib- 
rium for an initial state p7 

Theorem 2. Suppose that a Markovian system coupled 
to a heat bath is described by Lindbladian 7 with F- 

operators that form a set of generators for a nontrivial n- 
dimensional representation of a semisimple Lie algebra. 
The system returns to thermal equilibrium if and only if 
the representation is irreducible. 

Proof. Return to equilibrium for arbitrary initial state p 
will happen iff 

lim e^V = pa for all p e V{n), (8) 

where L = —i/h[HA, p] + C{p) denotes the right side 
of and Vili.) denotes the space of density operators 
on Ti. By assumption, the i^-operators are Hermitian, 
in which case ((S)) is equivalent to the statement that 
any operator X commuting with all the Va must look 
like X = c 1 for some constant c. This is equivalent by 
Schur's lemma to the statement that the Va generate an 
irreducible representation. □ 



D. Models of Decoherence 

Lidar et al. classified the subspaces in which 
there is generically no decoherence in this model, where 
"generic" means that the decoherence is suppressed 
independently of the noise parameters (i.e. the a- 
matrix of |5J)) and of the initial conditions. If g is 
semisimple (this assumption includes the Lie algebras 
of SU{n),SO{n),Sp{2n),G2,F4„Efi,ET,Es, and all di- 
rect sums of such algebras), then Lidar et al. P| char- 
acterize the decoherence-free (DF) subspaces as those 
which are annihilated by every one of the representa- 
tion matrices. At an intuitive level, this is not sur- 
prising since on such a subspace, the Lindblad operator 
^Fa,Ff) (p) — [FcpFjj] + [Fap, F^] is identically zero, and 
hence time-evolution is governed by the von Neumann 
equation. As representations of g, DF subspaces are al- 
ways composed of singlets. 

As discussed above, there are many interesting quan- 
tum systems which do not possess decoherence-free sub- 
spaces. Therefore, it is of interest to search for other 
effects which suppress decoherence to complement the 
existing techniques. On the other hand, at the highest 
level of generality, completely arbitrary decoherence is 
allowed by quantum mechanics. Thus it seems necessary 
to exploit symmetry in some way in order to achieve a 
suppression. Therefore we search for Lie-algebraic condi- 
tions which suppress the Lindbladian, but allowing some 
dependence on the noise parameters and on initial con- 
ditions. 

The noise parameters and initial conditions, while difh- 
cult to experimentally fine-tune, are also not entirely out 
of the hands of laboratory control. Further, because the 
solution operator e^* depends smoothly on any linear pa- 
rameters contained in L, the Lindbladian C{p) depends 
smoothly on the parameters aa.p, and hence a tiny vari- 
ation in aa,/3 results in a correspondingly small change in 
the dynamics. 

II. ALGEBRAIC IDENTITIES FOR LINDBLAD 
OPERATORS 

A. Suppression of Decoherence Without Noiseless 
Subsystems 

Before continuing, we record a set of simple algebraic 
identities for Lindblad operators. The Lindblad operator 
may be studied generally in terms of its matrix- valued 
bilinear form 

LA,B^2ApB^ -{p,B^A}. (9) 

Suppose that A, B are both Hermitian. We assert that 

La,b - 2[A,B]p + 2A[p,B]-B[p,A]-[p,B]A (10) 
La,b + Lb,a^[A,[p,B]] + [B,[p,A]] (11) 
La,a = [AAp.A]]. (12) 
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hence if A commutes with p then La. a = 0. If [p, B] 
commutes with A then 

La^b = 2[A,B]p + A[p,B]-B[p,A]. (13) 

The latter formula is antisymmetric in A, B. If [p, A] also 
commutes with B then (|13|l holds with A and B switched, 
hence La,b = ^Lb,a- This also follows immediately 
from (inj- If A] = A and [p, B] = B then 

La,b = 2[A,B]{p + l), hence La,b = -Lb,a- (14) 

We now investigate the consequences of Assumption 
n We identify elements a; G g with the operators 0(x) 
acting on TC. 

Lemma 1. Suppose that for each x G Q, [p,x] commutes 
with all of Q and Oap is real. Then C{p) = 0. 

Proof. The assumption entails that [p, Fa] commutes 
with Fjs for all a,p. Then Hill) implies that Lp^^p^ = 
—Lp^ p^ for all a, 13, but the sum L{p) involves these 
objects in symmetric combinations. □ 

It is natural to ask: when can the conditions of Lemma 
Qbe satisfied? If p can be expressed as a linear combina- 
tion of elements of g then the condition reduces to 

adp(0) C Z{q), 

where Z{q) is the center of q. 

Generally, it may not be possible to express p as a 
linear combination of (the identity and) generators of q 
[2^ . The vector space dimension dim(g) will often be 
smaller than — 1, the geometric dimension of the space 
of density operators. However, it is possible to equip q 
with an embedding into a larger algebra f), and note that 
the conditions of Lemma ^ will be satisfied if adp(g) is 
contained in the centralizcr of g within f), i.e. the set of 
all elements in (j that commute with all elements in q. 

Lemma 2. Suppose that adp is the identity on q (i.e. 
[yO, Fa] = Fa for all a) and aafj is real. Then C{p) — 0. 

Lemmas n and 121 point to the following general princi- 
ple: suppression of decoherence, for an initial state p, is 
related to the adjoint action of the density matrix p on 
the Lie algebra g containing the error generators. 



B. Lindblad Operators for su„ 

The above analysis leads to exactly calculable formulae 
for Lindblad operators in the n representation of su(n) . 
Let {Xj} be the canonical set of generators for su(n) in 
this representation, satisfying Ti^XiXj) = 2dij. Mac- 
Farlane et al. give a very clear discussion of su(n) 
generators and the identities that they satisfy; we use 
the same notation and conventions. 



The Lp^^Pi^ are sums of commutators, hence they are 
traceless and may therefore be expressed as a linear com- 
bination of the Xi. An arbitrary nxn density matrix p„ 
may be written in the Bloch representation 

a 

and write for Lxi,Xj- 

For calculating with Lie algebras, we find it extremely 
useful to use a convention even lazier than the Einstein 
summation convention to systematically not write sums 
and repeated indices, as follows. Each Xa is a matrix, but 
we never write the matrix indices; Va is a vector with the 
same dimension so v ■ X = VaXa is a matrix. Extend 
this to all tensors, so that v ■ fi, ■ X — > J2a c'^afabcXc 
etc. 

The standard basis for the su(7i) algebra satisfies the 
basic identities, which are derived in |25j: 

2i 2 

[PvjXb] = —V ■ fb- X, XaXb = -Sab^ + Qab ' X, 

n n 

[Xa,Xb]=2lfab-X, 

This makes explicit calculations easy. For example, if the 
terms in the diagonal standard form Q) happen to be the 
5u(n) generators, then using (|12|) we find 

2i 4 
La = [Xi, —v-fi-X]^-v-fi-fi-X. 
n n 

Using ,25|, eq. (2.12)], we then have 

J2Lti = -4.{vX)^A{l-np). (15) 

i 

Thus we have derived a very simple formula for the Lind- 
bladian of su(n)-decoherence. It vanishes iff p = ^1. 

III. SYMMETRY BREAKING 

Symmetry breaking in this context can be modeled as 
the introduction of new error generators which do not 
belong to the representation of g which gave the original 
symmetry. The standard assumption is that the sym- 
metry is broken perturbatively by modification of the 
system-bath Hamiltonian: 

Hab ^ Hab + eff/,sb . (16) 

where Hj sb denotes a symmetry-breaking interaction be- 
tween the system and the bath, and e <C 1. As discussed 
by Lidar et al. |^, the new terms added to the Lindbla- 
dian are of the form 

^'(P) = (""P-^^'ceGp + h.C^ + YbpqL,G^,eG, (17) 

a,p p.q 

where h.c. denotes the hermitian conjugate. Clearly the 
terms proportional to L^q^ ,^q^ are O(e^) and can be ne- 
glected in first-order perturbation theory, so we turn to 
the analysis of the 0{e) terms in the present context. 
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Definition 1. Let \) be a Lie algebra with subspaces g,f. 

(a) Following standard textbooks, we define the central- 
izer of g in t) by 

Z^(0) = {xei): [x,g] = V g e fl}. 

(b) Let : () ^ End(y) be a representation of I) on a 
vector space V and let p be an operator on V. We say 
that a Lie algebra homomorphism A : f) — s- f) switches 
centralizers for g, f if 

A{f)cZt,{g) and A(0)cZ„(f). 

We also say that p switches centralizers if A = adp does} 

Assume that f) admits a faithful representation (/):()—!■ 
g[(7i) on the system Hilbert space. Since (j) is faithful, 
identify f) with its image under 0, which makes sense 
since we are just interested in operators on Ti. 

Let f and q be the Lie algebras generated by the F- 
operators and the G-operators respectively. If p switches 
centralizers for f and g, then implies 

and therefore, the corresponding term in H17|) is propor- 
tional to 

Oap - a*ap = 2i lm{aap). 

If the couplings real, then the 0(e) symmetry 

breaking terms H17|) are zero, and the symmetry is per- 
turbatively stable (to first order). We have proven: 

Theorem 3. Let f and g be the Lie algebras generated 
by the F-operators and the G-operators respectively. If 
p switches centralizers for f and q, and Hap S R for all 
a,p, then C'{p) = 0. Ifaa.p ho-s an imaginary part, then 
the 0(e) terms in C'{p) are proportional to Im(aQ_p). 

Physically, these interactions are safe for quantum 
computing, because they do not contribute to decoher- 
ence. It is therefore of interest to locate examples of den- 
sity matrices which switch centralizers for certain subal- 
gebras. 

The most obvious example is the following. Suppose 
that f) = f g is a direct sum, so in an appropriate basis 
we can write all of the representation matrices in Block 
diagonal form; schematically. 




We will now continue with our convention of identify- 
ing f with its image under the faithful representation (j). 
Suppose that p is also Block diagonal in the same ba- 
sis, so that [p,f] C f and [p, g] C g. Under these con- 
ditions, the identity homomorphism switches centraliz- 
ers, i.e. f C .^i)(g) and g C 2'f,(f). Therefore, adp also 



switches centralizers. This, like any example with a DF 
subsystem, is a reducible representation, but it is plausi- 
ble that centralizer-switching occurs for irreducible rep- 
resentations as well. 

It is fruitful to apply this analysis to the case of a 
symmetry-breaking perturbation H17|l . Theorem |31 im- 
plies that if f and g switch centralizers, then the cross- 
terms, which are those proportional to eLp^.Q^ in H17|l . 
are zero. These are also the only 0(e) terms, so this kind 
of symmetry breaking is a purely O(e^) effect! 

IV. THE EFFECT OF COARSE-GRAINING 
A. Lie Algebra Channels 

Quantum channels defined by Lie algebra representa- 
tions were defined and studied extensively in a separate 
paper Here we recall the basic definitions and set 

notation, before discussing the relationship to Markovian 
dynamics. 

Definition 2. Let q denote a Lie algebra of dimension k, 
and let {Xi : i — I, . . . , k} be a basis of Q. Let (j) be a uni- 
tary representation of q on a Hilbert space Ti.. The gener- 
alized depolarizing channel or Lie algebra channel is de- 
fined to be the channel in which an error occurs condition- 
ally with probability p, causing an initial state {ip) Cz Ti. 
to evolve into an ensemble of the k .states 4'{Xi) , all 
with equal likelihood. 

The Kraus operators for the channel of Definition|21are 
given by 



Mo = 1, M, = 




where Z is a normalization constant, and will be deter- 
mined momentarily. The operators Mp_ are hermitian if 
the representation is unitary and if p G [0, 1]. 

Let K be the Killing form. If (j) is irreducible and the 
chosen basis satisfies pseudo-orthonormality 

KiX„X,) = n6^,, n^O, (19) 

then J2z(l^i^i)'^ = Z ■ 1 with Z = nc2(0) where C2((/)) 
is the quadratic Casimir, and therefore the condition of 
probability conservation MpM'^ = 1 is satisfied. If (j) is 
irreducible but the basis does not satisfy (|19l) . then Dcf.|21 
does not define a channel. If (f) is reducible, then one must 
replace Z Zp where the Zp are the normalization 

constants for the independent irreducible components. 

Using H18|) . the Lie algebra channel has the explicit 
Kraus decomposition 

k 

p^£{p)^{l-p)p+^Y. '^(^') P '^(^') ■ (20) 



^ While Def. 0|a) is standard, I have never seen Def. 0^b) before. 



In what follows, we will write Xi for (t>{Xi) since the 
distinction is clear from context. 
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The companion paper [23| gave methods for calculat- 
ing with these channels using Lie algebra identities. Let 
Eg denote the channel based on the fundamental repre- 
sentation of g; for example, i£'su(n) is based on the n of 
su{n). The action of £ciu{n) on a general hermitian ma- 
trix p — n^^(tr(p)l + V ■ X) can be calculated exactly in 
closed form, yielding 

£{p,) ^ - f tr(p)l + i^-P)^'-^ y.x). (21) 
n \ — 1 y 

In the qubit case, the coefficient of w is 1 — 4p/3, which is 
exactly consistent with the standard qubit depolarizing 
channel ©-Q. 

B. Coarse-graining 

Markovian dynamics may be derived from quantum 
theory of measurement by coarse-graining time. Ex- 
panding Lindblad's equation p — — i aa{[Fa, pF^] + 
[FaP, F^]) to first order in the short time interval r yields 

p{t + r) = MQp{t)Mo + T ^ Fa,p{t)Fi , 

with 

Suppose that the Lindblad operators equal to properly 
normalized Lie algebra generators, Fi = Z^^I'^Xi. Then 

A/o = (l - 1 ~ 1 + 0(r2). 

Also defining 




shows that the coarse-grained Lindblad's equation is 
completely equivalent to the Lie algebra channel H18|) to 
first order in t, under the mapping t —* -p. 

It is interesting that in the Lie algebra channel, p is 
probability of error while in Markovian dynamics, r is 



the coarse-graining time. In fact, their equivalence to 
first order coincides with one of the defining properties 
of a Poisson process: that the probability of exactly one 
change in a sufficiently small time interval of length r is 
proportional to r, with constant of proportionality given 
by the event rate for the generation of errors. This as- 
sumption is often made (either implicitly or explicitly) 
in the theory of quantum error correction and describes 
a quantum Poisson process. An example of a quantum 
Poisson process is spontaneous emission |2^. For such 
processes, we expect p oc t to first order in p or r. 

If we specialize to g = su(n), then the evolution of 
the channel was found to be (PT|l with p replaced by r. 
Since H21|) is already linear in p, any O(t^) terms just 
come along for the ride. The conclusion is that in su{n) 
Markovian dynamics, 

trp(t) 1 / rn^ \ , ^ , ^ 

n n\ — 1/ 

Eqn. I|22|l is something which, a priori we did not expect: 
an explicit formula, valid to first order in the coarse- 
graining time, for the nontrivial Markovian system de- 
fined by 5U„. The simplification here is related to the 
simplification (|15(l of the Lindbladian. 



V. CONCLUSIONS 

Since there are many examples in which no 
decoherence-free subsystems exist (among them all cases 
where the error generators act irreducibly on the system 
Hilbert space) , it is of interest to search for novel mecha- 
nisms which suppress decoherence in these more general 
cases. We presented three results (Lemmas^andEl The- 
orem which indicate decoherence suppression without 
the need for noiseless subsystems. There is a certain 
trade-off; our methods do not necessarily apply to an ar- 
bitrary initial density matrix, or for completely generic 
noise parameters aaf^. On the other hand, our compu- 
tational methods are novel and the result-suppression of 
decoherence in the error-algebra approach without noise- 
less subsystems-is an interesting new direction which 
warrants further study. 
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